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Abstract 

Two-point space-time correlations of air density fluctuations in unheated, fully-expanded free jets at 
Mach numbers Mj = 0.95, 1.4, and 1.8 were measured using a Rayleigh scattering based diagnostic 
technique. The molecular scattered light from two small probe volumes of 1.03 mm length was measured 
for a completely non-intrusive means of determining the turbulent density fluctuations. The time series of 
density fluctuations were analyzed to estimate the integral length scale L in a moving frame of reference 
and the convective Mach number M c at different narrow Strouhal frequency (St) bands. It was observed 
that M c and the normalized moving frame length scale L*St/D, where D is the jet diameter, increased with 
Strouhal frequency before leveling off at the highest resolved frequency. Significant differences were 
observed between data obtained from the lip shear layer and the centerline of the jet. The wave number — 
frequency transform of the correlation data demonstrated progressive increase in the radiative part of 
turbulence fluctuations with increasing jet Mach number. 

I. Introduction 

The statistical description of turbulence is usually based on the two point space time correlations of 
flow fluctuations (ref. 1). Various correlations commonly employed in turbulence modeling are subsets of 
information obtained from the general space-time correlation measurements. Therefore it is natural to see 
that the description of the noise sources in the acoustic analogy (ref. 2) or other schemes (ref. 3) depends 
on the space-time correlation of the turbulent field. All existing jet noise prediction schemes, that try to 
model the underlying physics, bank on experimental data to describe the spatial and the temporal 
variations of the correlation coefficients. The need for the correlation data was realized as soon as the 
LighthilTs celebrated work was published and hot-wire measurements by Laurence (ref. 4), Davies et al. 
(ref. 5) and others provided insight. Like the prior modeling efforts most of these experimental data treats 
turbulence as all random fluctuations with single length and time-scale describing fluctuations at all 
frequencies. It has been realized that further improvements in noise prediction can not be achieved 
without considering the differences at different frequency ranges. The large-scale/fine-scale approach of 
Tam et al. (ref. 6) and the frequency dependent scales approach of Harper-Bourne (ref. 7) are examples of 
this change in the modeling approach. The present experimental effort is motivated by this new approach. 
The goal is to determine frequency dependence of parameters such as convective velocity, integral length 
scales, and time scales that are measured via space-time correlation. 

The first significant measurement of space-time correlation in free jets used hot wire probes (refs. 4 
and 5); the jets under study were in the Mach number range of 0.2 to 0.55. To avoid excessive 
interference from probe supports inserted into the flow, usually measurements were done only in the 
peripheral shear layer. Nevertheless, the observations made in these tests are still used in noise source 
modeling. Some of the significant observations are the following, (a) The integral length scale along the 
axial direction in a fixed frame of reference varies linearly with the downstream distance ( x = 0.1 3x for 

x/D<6. (b) The radial length scale is nearly a third of the axial one: ( r = ( x /3. (c) The radial variation of 

£ x is small, (d) The convective velocity Uc varies from 0.7Uj at the inside edge of the shear layer to 0.3Uj 
at the outside edge with around 0.65Uj at the max shear location. Uj is the jet centerline velocity, (e) The 
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integral time scale in a moving frame of reference, calculated from the time envelope of correlation with 
different probe separation, is inversely proportional to the local shear: T x ocl/(5U/5r), where U is the local 
axial velocity and r is the radial direction. The integral length and time scales are defined later in the text. 
An interesting debate started by these early works in 1960s’ was whether to use an exponential decay 
function or a Gaussian decay function to model the time variation of correlation. It has remained an issue 
even today (refs. 8 to 10). Since all of these observations were made in low speed jets and primarily in the 
lip shear layer, there is a need to determine variations in other regions of the jet and also in the practical 
Mach and Reynolds number range. For example, according to the observation (e), moving frame integral 
time scale becomes unphysical at the jet centerline where shear gradient is zero. 

Early on Fisher and Davies (ref. 11) recognized the need for determining the frequency dependence 
of the turbulent scales. They band-pass filtered hot-wire signals before performing correlation. A 
drawback of this approach is that the filter’s impulse response causes spurious oscillations in the filtered 
correlation data. The true temporal behavior of the correlations, especially at longer time delay, becomes 
difficult to isolate. The narrower the frequency band the more persistent is the oscillations. The present 
work takes a frequency domain approach used by Harper-Bourne (ref. 7) that is more conducive to 
investigating frequency dependent correlation. 

The modem optical techniques such as LDV (refs. 12 and 13), PIV (ref. 9) and optical deflectometry 
(ref. 14) overcome many limitations of hot-wire probes. The advanced computations based on LES 
(ref. 15) and DNS (ref. 16) methodologies have also joined the force to determine turbulence scales, but 
their efforts are mostly limited to low Reynolds number flows and over limited frequency ranges. The 
present Rayleigh scattering based diagnostics allowed for a completely non-intrusive means of measuring 
turbulent fluctuations. There are no seed particles used in the technique, instead light scattered by the gas 
molecules are analyzed. In principle, the scattered light can be used to determine velocity, density and 
temperature fluctuations, however; only density fluctuations are reported in the present paper. Density 
and axial velocity fluctuations are found to be well correlated; hence the space-time correlations measured 
from the density fluctuations are expected to be similar to that of the axial component of velocity. 
Additionally, noise sources can be described via any flow variable and for the supersonic jets density is 
easily measurable. 

The present work is a part of a large database development project for validation of the CFD and 
CAA codes. The time averaged velocity, density and temperature profiles, density fluctuations spectra, far 
field noise spectra, and noise sources identified by direct correlations between far field noise and flow 
fluctuations were reported earlier in references 17, 18, and 19. 

Basics of Space-Time Correlation 

The two point space-time correlation is defined either in a fixed frame of reference or in a moving 
frame that is convected with the flow at velocity U c . When a probe is kept at a fixed position of the jet the 
measured temporal fluctuations are in a fixed frame of reference. The direct correlation measurement 
from two such probes placed at X and X+r| (q is the separation) provides fixed frame correlation. Note 
that in general, the probe separation can be in the axial, radial or azimuthal direction, however; for the 

following discussion separations only along the axial direction is assumed so that q = i q x where i is the 

unit vector in the stream wise direction and q x is the stream wise component of q. For the present case of 
density fluctuation it is defined as the following: 

1 T 

(p' * P') (X, q x , T) = — J p'(X, t) p'(X+ iq x , t +x) dt (1 ) 

1 0 

Where, p'(t) is the fluctuating part of the instantaneous air density obtained by subtracting the average: 
p'(t) = p(t) - p, and x is the correlation time. If the probes were allowed to move with the local convective 
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velocity, the resulting correlations would have been in a moving frame. Since turbulent fluctuations are 
convected by the flow, the signal from a fixed probe is expected to show rapid fluctuations and is 
expected to de-correlate over smaller spatial and temporal separations: q x and x. However, if the effect of 
the convection is avoided by measurements in a moving frame, then correlation values represent the 
decay/growth and the distortion of eddies. This is expected to occur over longer separations. In spite of 
the modeling advantages (ref. 2), it is difficult to perform correlation measurements in a moving frame. 
Therefore, the moving frame properties are approximated from the fixed frame measurements. Note that 
in a moving frame the space-time correlations are related by convective velocity. Hence, if the experiment 
is conducted in a laboratory fixed frame the moving frame auto-correlation can be approximated by 
selectively choosing peaks in fixed-frame correlation (ref. 5): 


(p'*p') c ( x ,x) = (p'*p')( x ,r| x =U c x,x) (2) 

The superscript c denotes the moving frame correlation. The correlation values are normalized by the 
standard deviation of fluctuations p rms at the two measurement points and expressed as correlation 
coefficients: 


r(X,r| x ,T)= 


(p'*p')(X,r| x ,x) 

Prms(X)O) P rms (X, %)' 


(3) 


The definition of the convective velocity has been a subject of some discussion in the past 
(Wills (ref. 20)). It is defined as the ratio of r|/x that maximizes the correlation coefficient r. The 
derivative can be taken with respect to either r) or x. The latter is usually the case for a fixed separation r| f ; 
we look for x c such that: 


Maufi 0, U c =2t (4) 

fix X c 

It also can be defined in terms of a fixed time delay Xf that maximizes the correlation (Willmarth and 
Wooldridge (ref. 21)): fir(r|,Xf)/<3r| = 0 and U c = r| c /x f . Since in both situations, the convection velocity was 
found to be somewhat dependent upon the choice of the fixed parameter Wills (ref. 20) introduced a 
global average. The present work follows the commonly used definition of equation 4 to calculate the 
convective velocity. 

For a generalized model of space-time correlation that can be incorporated in the turbulence (noise 
source) modeling, there is a need to determine length and time scale of turbulent eddy. This is done via 
integral length and time scales (Townsend (ref. 1)). In a fixed frame of reference, the integral length scale 
t and time scale t T are defined as: 


4x)=|r(X,T lx ,0)dT 1 x , t T (x)=Jr(X,0,x)dx (5) 

0 0 

As discussed earlier integral scales defined in a moving frame (length L, time T x ) provide better insights 
into the distortion of turbulence. This can be calculated from the moving frame auto-correlation function 
as follows: 
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(6) 


T x (x)= | r c (X,x)dx =Jr(X,ri x = U c x,x)dx, L(x) = T x U c 
0 0 

For a frozen pattern of turbulence, r c = 1 and L = go. However in realistic flows r c — > 0 as r | x — > <» and both 
L and T x are finite. As discussed earlier the moving frame scales are longer than those in a fixed frame, 

T x » t x . 

To introduce frequency dependent scales, the correlation function in the equation (1) is Fourier 
transformed in time to obtain cross spectrum G xr| : 


GO 

G XT1 (X,q x ,co)= |(p'*p')(X,q x ,x)e i ^dx = G x G; = |G XT1 |e- i<I> (7) 

—00 

Here G x and G n represent Fourier transform of density fluctuations at X and X+iq x , and the superscript * 
represents complex conjugate. The cross-spectral phase ®(co) reflects the time shift x c in fluctuations 
between the two measurement points and can be expressed as: ® = co x c . The primary cause of the time 
shift is the convection of turbulence fluctuations by the mean flow x c = r| x /U c . Any distortion of the 
turbulent eddies over the separation distance is also manifested in the phase. However, for small probe 
separation the latter contribution is expected to be small. Therefore, the frequency dependent convection 
velocity can be estimated as: 


U c (x,eo) 


Tx® 

® 


( 8 ) 


To reduce the measurement uncertainty, and to somewhat cancel out the contribution from the distortion 
of eddies, the convective velocity is calculated via averaging over multiple separation points. 

Returning back to equation (7), after stripping away the phase part (and thereby the effect of 
convection), the remaining magnitude part of the cross spectrum is indicative of the change in correlation 
seen by the probes if one of them were convected with velocity U c . This approximates the decay in 
correlation in a moving frame of reference. Hence the magnitude of the cross-spectrum can be used to 
estimate the integral length scale in a moving frame of reference. Towards this goal, the cross-spectrum 
needs to be normalized using the individual auto-spectra, and a coherence function is introduced: 


R(x,r| x ,co) = 




G X r| 

e 

,2V /2 A 

1? 1 

_ ,2V /2 

o 

x| G n 


where, |G X | 2 = G X G*, |G in | 2 =G n G* 


( 9 ) 


Harper-Bourne (ref. 7) used a coherence function which is slightly different from the above expression. If 
|G X | 2 ~ |G n ; 2 or magnitudes of the turbulence fluctuations are similar then: 


R(X,r| x , co) ! 


J XT) 


|G, 


e - i0 =Re -io 


( 10 ) 


R(X, r| x , co) is a simplified coherence spectrum (described as coherence modulus by Harper-Bourne) and 
is indicative of frequency dependent correlation coefficient. In general R is unity for perfect correlation 
and zero for complete de-correlation. A downside of the simplification is that in a growing shear layer 
with continuous turbulence production R can take values greater than unity. This was frequently 
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encountered in the present data. The frequency dependent scales in a moving frame can now be defined 
as: 

L(x, co) = J R (X, q x , cd) dq x , T t (x,co) = ^&^ (11) 

o c 

These integral scales and other parameters defined above are measured in the present work. 

A great utility of the space-time correlation data is that it allows for a direct determination of the part 
of turbulent fluctuations capable of radiating to the far field. Sound generation process involves a 
coupling between the fluctuations measured in the source region and the ambient medium. Most of the 
fluctuations in the source create only evanescence waves that decay down within a few diameters from 
the jet. Only a part of the fluctuating energy that attains the speed of sound along an observer direction in 
the ambient medium is ultimately radiated to the far-field (ref. 22). A turbulent eddy convected at 
subsonic speed and without any distortion does not radiate sound. However, distortions associated with 
growth and decay enables a subsonic eddy to radiate (ref. 23). To determine this radiative component of 
fluctuations, one needs to perform a wave number (k) - frequency (co) transform of the correlation data. 
Towards this end at first correlation data (p' * p')(X, q x , x) are measured from a large number of 

separations. Subsequently, space-time Fourier transforms shows the distribution of energy in the k-co 
space: 


00 00 

S k< a(X,k, ©)= | J(p'*p')(X,q x ,T)e i “ T e ikl ixdTdq x 

—00 —00 


( 12 ) 


In the far field, away from the jet, k-co spectrum is a straight line co/|k| = c x , where c x is the speed of 
sound in the ambient medium. In contrast, due to the dispersive nature of turbulence, where fluctuations 
at a given frequency may have a range of convective velocity and therefore many different wave numbers 
k, k-co spectrum inside the plume is significantly broadened. The part of the spectrum with co/|k|>c x is 
capable of coupling (radiating) to the far field at different observer angles. Note that the wave-number is a 
vector quantity with components along the axial, radial and tangential directions. In the present work the 
S kt0 spectra in axial wave-numbers were calculated for different jet Mach numbers and at different 
locations inside the plumes. 


II. Experimental Setup 

The experiments were performed in an unheated small jet facility at the NASA Glenn Research 
Center. A convergent nozzle was used to produce subsonic plumes and convergent-divergent nozzles 
were used for the Mach 1.4 and 1.8 plumes. All three nozzles had a nominal equivalent exit diameter of 
D = 2.54 cm. A detailed description of air density measurements using Rayleigh scattering principles can 
be found in reference 18. Fundamentally, Rayleigh scattering is the result of elastic scattering (no change 
in frequency between the incident and the scattered light, except for the Doppler shift) of light caused by 
molecules, such as the Oxygen and Nitrogen molecules, present in the air. The intensity of the scattered 
light is directly proportional to the molecular number density, and therefore the bulk density p at the 
probe location. For a fixed composition gas, such as air, and for a fixed optical setup the changes in the 
scattered light intensity is directly related to the changes in the air density at the probe location. This 
intensity measurement principle is used to determine turbulent density fluctuations in the present setup. 
Measurements using the Rayleigh scattering principle do not require any seed particles, and therefore, the 
nuisances and the uncertainties involved in particle based measurement techniques are avoided. On the 
downside, to avoid large scattering from the contaminated particles special precautions have to be taken 
to cleanse the jet of dust particles. Reference 18 provides a detail description of the steps taken towards 
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this goal. For the sake of brevity only the new part, created for the two beam setup, is described in the 
following. 

The present optical setup around the nozzle facility is different from the earlier version used to 
measure single-point statistics. Figure 1 shows schematics of the facility and figure 2 shows a photograph 
of the facility. The optical setup consists of two beams crossing the jet plume in vertical and horizontal 
planes. The top view in figure 1(a) highlights the optics for the horizontal beam while that of the vertical 
beam is better seen in the side view of figure 1(b). Laser light scattered from a small region of each beam 
was collected using two separate sets of lenses oriented 90° to the laser paths, and measured by photo- 
multiplier tubes. The scattering angle is of little importance for density measurement; however; the setup 
was built with an ambition of simultaneously measuring density and velocity fluctuations from two 
different points in the jet. Since velocity measurement involves determination of the Doppler shift of the 
scattered light, the beam path and the collection direction are designed such that the same axial 
component of velocity could be measured by either setup. 

Two mounting frames holding beam- forming optics, beam dump and collection optics for each beam 
were mounted on separate 2-axes traversing units. The traversing units allowed independent positioning 
of the probe volumes in an axial(x)-radial(r) plane. The light source is a continuous-wave, frequency- 
doubled Nd:V0 4 solid state laser operating at 532 nm wavelength. The laser head was placed at the 
bottom of the structural frame attached to the horizontal traversing unit. The primary beam was split into 
two parts by a 30 to 70 beam splitter. The weaker 30 percent part of the beam was passed directly via a 
focusing lens and a half-wave plate to form the horizontal beam. The horizontal probe volume was set at 
the focal point where the beam waist is about 0.16 mm. The beam is ultimately dissipated in a beam 
dump. The second, 70 percent part of the primary beam was coupled to a multimode, 100 micron core 
diameter optical fiber and passed to the top of the frame holding optics for the vertical probe volume. A 
combination of two lenses, 80 and 300 mm focal lengths, focused the light emerging from the fiber to a 
0.38 mm waist, which constituted the vertical probe volume. The collection optics for each unit were 
made of a combination of 160 and 300 mm focal lengths and 80 mm diameter lenses, which focused the 
collected scattered light from each probe volumes on a 0.55 mm core diameter multimode fiber. The fiber 
diameter fixed the probe volume length (extent of the laser beams imaged on the fiber) to 1.03 mm. 
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Figure 2. — Photograph of the facility. 

Attention was paid to minimize vibration on the optical components by the high noise level produced by 
the jet. The laser head was enclosed in a semi-anechoic chamber and mounting frames were placed 
upstream while the laser beam propagated downstream. 

Since the intensity of the Rayleigh scattered light is polarization dependent, the horizontal beam was 
vertically polarized to maximize scattering towards the horizontal collection optics. The vertical beam 
was passed through a multi-mode fiber which scrambled the incident polarization, and therefore; 
scattering from half of the vertical beam was effectively collected by the vertical collection optics. The 30 
to 70 splitting of incident laser power was guided by the desire to achieve nearly the same scattering 
power at the two collection setups. Indeed at the end, when the probe volumes were separated, similar 
rates of photo- electron counts were achieved. Flowever, when the beams were made to overlap the 
differences in polarizations caused the photon arrival rate to the horizontal collection optics to exceed that 
towards the vertical. 

The scattered light, collected by the optical fibers, was passed to an adjoining room where two 
photomultiplier tubes and photo-electron counting electronics were kept (not shown in figs. 1 and 2). 
According to the Rayleigh scattering principle, for a fixed molecular composition gas (air in this instance), 
the intensity of the scattered light is directly proportional to the local air density. Therefore, the photo- 
electron arrival rate is related to the air density at the probe volume. The proportionality constant was 
determined through an in-situ calibration (ref. 1 8). It was performed in the plume of the convergent nozzle 
operated in the subsonic Mach numbers. At each operating condition the photon arrival rate was counted and 
the jet density was calculated using isentropic relations. Subsequently, a straight line was fit through the data 
to determine the proportionality constants a and b. 
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N = (a p + b) At 


(13) 


The additional constant b is needed to account for the room light and stray scattered laser light. Since two 
counters were used, two sets of calibration constants ai, bi, and a 2 , b 2 were calculated. 

As described earlier, the cross correlation data are normalized by the auto-correlation at the reference 
location where both beams overlap. To avoid the electronic shot noise contribution (this will be described 
later in the text) auto-spectra are calculated by cross-correlating counts from two PMTs with overlapped 
probe volumes. Since the effective laser power nearly doubles when the beams overlap, a new set of 
calibration constants were needed to convert the photo-electron arrival rate to physical density numbers. 
Therefore, a total of four sets of calibration constant were calculated: two sets (one for each beam) with 
separated probe volumes and two with overlapped probe volumes. Overlapping the probe volumes involved 
maximizing the count rates by trial and error. This was achieved relatively easily at one point in the plume. 
However; due to a slight misalignment of the two traversing units, a problem appeared in maintaining the 
perfection at all other locations. The misalignment was small: few tenths of a millimeter over the traversing 
range, yet sufficient to cause a drift in calibration. Therefore, the overlapped probe calibration was adjusted 
based on the known density (measured by one of the separated probes) at the overlapped point. 

The independent x-r traversing units allowed arbitrary separation of the probe volumes. The minimum 
separation, that avoided “cross-talk” between the measurement points, was determined by the diameter 
(0.55 mm) of the optical fibers onto which the beams were imaged, the scattering angle (90°), and the probe 
volume diameters (0. 16 and 0.38 mm). The laser beam profile, at the focal point, is Gaussian and the 
diameters simply refer to the full width at half maxima of this profile. The tail of the laser beam extends 
beyond the prescribed diameter. Hence, it was found that the minimum separation for which independent 
measurements by the two probe volumes can be achieved is about 2 mm (0.08D). 

Returning back to the signal processing, to resolve time variation in density, photo-electron counting was 
performed in a series of contiguous time bins typically of At = 8 192 ns duration (the bin duration was 
changed for different Mach number jets). This produced an effective sampling rate of 122 K samples/sec. The 
photon counting process also introduces an effective low-pass filter whose description can be found in 
reference 18. Usually two long records, each n tot = 5,120,000 points long, were collected for every 
measurement. The standard deviation of density fluctuations at each probe location was estimated from the 
power spectrum of the photon counts. The average values need to be subtracted from the individual strings, 
N/= Ni - N av (i = 0, 1, 2, . . .. n-1), before the power spectrum is calculated: 


P N / 2 (f 1 ) = 2 F N /(l) , f 1 = 


n-l 


nAt 


1 = 0,1, 2, — -1 where, F N /(1) = — exp 


n 


,27til' 


(14) 


i=0 


Here, f is the frequency associated with each line. The Welch method of modified Periodograms (ref. 25) 
was used to calculate the spectral density. Each long record was divided into small segments of n = 5 12 data 
points. The adjacent segments were overlapped by 50 percent. The modified periodograms of corresponding 
segments from the two PMTs were calculated and then used to determine local estimates of spectral density. 
All local estimates were averaged to obtain the final spectral density. The Power spectrum of the air density 
fluctuations requires multiplication by calibration constants: 

P p” <f,) = ^ P N» <f,) ( 15 ) 


The problem with this straightforward method is that the calculated spectrum becomes overwhelmed from 
the electronic shot noise, an unavoidable consequence of low light level measurements. Shot noise poses the 
fundamental uncertainty in all light intensity measurements. Following Parseval’s identity: 
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1 = 0 



+ CT 


2 

P- 


In other words, the mean-square calculated from the spectrum is a sum of shot noise a^ h and desired air 

density fluctuations . An estimate of the shot noise contribution can be obtained by noting that it is white 

noise and contributes equally to all frequency bins. The constant shot noise floor was estimated as the 
average of spectral density values at the highest end of the resolved frequency: of 55 to 61 KHz 
(Nyquist limit): 


Sshot : 


1 

'6000 


61000 

J P N / 

f =55000 


(17) 


The energy from turbulent fluctuations is expected to be small over this frequency range. 
Since a^ h = [ S s [ 10 t df , the standard deviation of density fluctuations is: 


Prms 



^shot ] 


(18) 


Individual power spectral density and shot-noise floors in PMT spectra are shown in figure 3(a). The data 
presented are in power spectral density (rather than power spectra) where the power in the individual bins are 
divided by the line width 


Gl(X,f,) 


AtP p Q 

2n 


(19) 


and frequency is normalized to Strouhal number St = fD/Uj; Uj is the velocity in the potential core. The cross- 
spectral calculation followed a similar path where the Fourier transform of one data series was multiplied by 
the complex conjugate (indicated by 3 * in the following equation) of the other: 


G xll (X,q x ,f|) = 


!>[(■> F t;d>) 

&l &2 Atn 


( 20 ) 


The magnitude of the cross-spectral density is shown in figure 3(a) and phase in figure 3(c). The frequency f 
defined in this equation is related to the circular frequency in equation (7) via co = 2nf. The cross-spectral 
function, for the discretely sampled data, defined in the above equation is the stepping stone towards 
calculating correlation parameters defined in the Introduction. For example, the magnitudes of cross-spectra 
from a set of data with different separation q x is used to calculate frequency dependent integral length scales 
while the phase variation is used for determining the convective velocity. 

Notice that when the probe volumes are coinciding, i.e., q= (0, 0, 0) the cross-spectrum becomes the 
auto-spectrum of the density fluctuations at the probe location X. Interestingly this is the true auto-spectrum 
without the influence of electronic shot noise. Shot noise from two photo-multiplier tubes are uncorrelated; 
hence a cross-correlation between the two signals rejects this uncommon part, bringing out the desired 
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Figure 3. — (a) Auto-spectrum of individual PMT signal and magnitude of density cross-spectrum 
(b) cross-correlation and (c) cross-spectral phase in Mj = 1.4 jet with one probe at (3.0D, 0.5D) 
and the other at (3.5D, 0.5D). 

spectrum of air density fluctuations. In reality, due to a finite averaging time, a noise floor remains as seen in 
the cross-spectrum presented in figure 3(a). This floor ultimately limits the usable frequency range for the 
correlation study. The cross-correlation was calculated via inverse Fourier transform F _1 of the cross- 
spectrum: 


(pi*P2)(x,t1x,t) = - 


f N ( (1) • 


F Ni® 


ai a 2 


At 2 


1 = 0, + 1, ±2, ± — 1. 

2 


( 21 ) 


Figure 3(b) shows one such cross correlation plot. Following is a discussion of the uncertainty in the 
normalized correlation coefficients. 

The cross-correlation process mostly eliminated contribution from the electronic shot noise. Additionally 
the large data sets used for calculation produced fairly converged statistics. Flowever, the correlation values 


are normalized by individual standard deviations, 


(pl*P2)(x,r|xA) 


which are the primary sources of 


Pi rms (X) P 2 rms (X,Px) 

uncertainty in the presented form. Calculation of the standard deviation of density fluctuations required shot 
noise subtraction as described earlier. The subtraction process produces estimation of the standard deviation 
within ±5 percent uncertainty. Flence, a maximum uncertainty level of ±10 percent is expected on the quoted 
normalized correlation values. This large uncertainty supersedes all other sources, such as the contamination 
from stray particles, numerical convergence error etc. It is also expected that the highest uncertainties are 
present in the lowest Mach number condition where the absolute value of density fluctuation is the lowest. 
The uncertainty level decreases with an increase in the jet Mach number. 


Results and Discussion 

Profiles of time averaged density and standard deviation of density fluctuations for the same jets under 
study were reported earlier in reference 1 8, and will not be repeated in this paper. The space-time correlation 
data were measured from many different locations in the plumes. Figure 4 shows distributions of the density 
fluctuations in the plumes and also identifies reference locations about which correlation data were obtained. 
The reference position signifies the location of the fixed probe volume. The second probe volume was mostly 
separated along the axial direction and, in a few cases, along the radial direction. Figure 4 demonstrates that 
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Figure 4. — Distribution of normalized density fluctuations in the indicated Mach number jets. 
The white dots indicate reference positions about which space-time correlations were 
measured. 


an increase in the jet Mach number causes a progressive reduction in the growth rate of the lip shear layer 
which manifests in a lengthening of the potential core. Note that jets under study are unheated, i.e., the total 
temperature of the plume is the same as that of the ambient. The change in density is a manifestation of the 
compressibility effect. Therefore, the density fluctuations discussed in this paper are a reflection of the 
fluctuating compressibility which is expected to be closely related to noise sources. 

The bulk of the figures presented in this paper are from two sets of reference points. For the first set, the 
reference point was kept fixed at the lip line and three diameters from the nozzle exit X(3.0D, 0.5D). For the 
second set the reference point was placed at the centerline and three different axial positions x/D = 8, 10, and 
12 for the three plumes of Mj = 0.95, 1.4, and 1.8, respectively. The difference in the axial location is a 
reflection of the increase in the potential core length Lp, and nominally represents an equal distance of 
~1.2 Lp. 

Figures 5(a) and (c) shows two sets of normalized correlation plots for the reference probe located at the 
lip shear layer and at the centerline. The temporal dependence of the correlation curves are of interest for 
noise source modeling. Interestingly, data from the shear layer and the centerline show marked differences; 
the former shows significant oscillation which is absent in the latter. Many jet noise prediction schemes 
(refs. 8 and 10) model the temporal dependence by either an exponential decay: exp (-x/io) or a Gaussian 
decay: exp (-x 2 /x 2 0 ) function. Clearly the centerline data shows better compliance to this modeling assumption 
than the shear layer counterpart. Additionally a closer scrutiny of the data shows that profiles become more 
asymmetric as the separation between the probe locations is increased. 

To calculate the integral length scale in a fixed frame of reference (eq. 5) correlation coefficients at 
x = 0 need to be plotted as a function of the probe separation q x . Figures 5(b) and (d) show two such plots 
from two reference locations. Determination of the integral length scale involves calculating the area under 
the r-q curve. While that is feasible for the centerline data, the negative oscillation in lip line data of 
figure 4(b) is expected to produce very small integral length scale. A length scale is a measure of the spatial 
extent of an eddy. In the lip shear layer, correlations persist over longer separations, indicating longer length 
scales; not the small value expected from figure 4(b). Therefore, calculations of the integral length scale in a 
fixed frame of reference were not pursued. This problem is avoided in calculations made in a moving frame 
of reference via the use of coherence spectra. The advantages of the frequency dependent scales were used in 
the recent noise modeling efforts of Goldstein and Khavran (ref. 24). 
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Figure 5. — (a), (c) Time domain correlation with the indicated probe separation 
and (b), (d) variation of the correlation coefficients at zero time separation 
from two different reference locations in Mj = 0.95 jet. The reference probe is 
at (a), (b) lip line X(3D, 0.5D) and (c), (d) centerline X(8D, O.D). 


Envelopes of the correlation curves in figures 5(a) and (c) show the decay in the turbulent density 
fluctuation observed in a moving frame of reference. The integral time scale (in the moving frame) can be 
calculated from this envelope, representing the auto-correlation function in a moving frame, using 
equation (6). An exponential decay function was used to fit each envelope and the integral time scale T x was 
calculated. The normalized time scales are presented in table I as “wideband” values. The time domain 
correlation data of figure 4 does not explicitly show frequency dependence of integral scales; T t provides an 
average property over all frequencies. The time separation between the peak values of the correlation is a 
measure of the average convection time for an eddy to move from one probe location to the other. This 
allows for an estimation of the convective velocity, which once again is an average of all frequencies of 
fluctuation. Such wideband convection velocities are shown in table II. After knowing U c the wideband 
integral length scale in a moving frame can be calculated from the integral time scale: L x wide /D = T t U c /D. 
This is also shown in table I. 
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TABLE I.— INTEGRAL SCALES IN MOVING FRAME OF REFERENCE 


Mj 

Reference probe 
location 

Time scale 
T t Uj/D 

Integral length scale L x /D 


x/D 

r/D 

Wideband 

Wideband 

St = 0.3 

0.7 

1.0 

1.6 

2.5 

3.5 

0.95 

3.0 

0.5 

1.29 

0.9 

2.88 

1.44 

0.67 

0.41 



0.95 

8.0 

0 

1.53 

1.18 

2.01 

0.99 

0.85 

0.65 



1.4 

3.0 

0.5 

1.22 

0.88 

3.01 

1.81 

1.16 

0.46 

0.33 


1.4 

10.0 

0 



1.84 

1.03 

0.84 

0.69 



1.8 

3.0 

0.5 



1.85 

2.12 

1.38 

0.7 

0.44 

0.33 

1.8 

12.0 

0 



2.2 

1.07 

0.85 

0.65 

0.48 

0.43 

0.95 

3.0 

0.4 



3.27 

2.26 

1.19 

0.76 



0.95 

3.0 

0.6 



2.15 

0.93 

0.45 





TABLE II— CONVECTIVE MACH NUMBER: M c = U c /Uj 


Mj 

Reference probe 
location 

Convective Mach no. M c 


x/D 

r/D 

Wideband 

St = 0.3 

0.7 

1.0 

1.6 

2.5 

3.5 

0.95 

3.0 

0.5 

0.7 

0.63 

0.72 

0.75 

0.74 



0.95 

8.0 

0 . 

0.77 

0.73 

0.78 

0.82 

0.88 



1.4 

3.0 

0.5 

0.72 

0.74 

0.77 

0.78 

0.76 

0.77 


1.4 

10.0 

0 


0.76 

0.81 

0.83 

0.88 



1.8 

3.0 

0.5 


0.66 

0.77 

0.8 

0.82 

0.86 

0.86 

1.8 

12.0 

0 


0.74 

0.8 

0.83 

0.88 

0.9 


0.95 

3.0 

0.4 


0.67 

0.75 

0.78 

0.76 



0.95 

3.0 

0.6 


0.58 

0.7 

0.62 





As described in the Introduction, the first step towards calculating frequency dependent scales is to 
determine auto- and cross- and coherence R (eq. 12) spectra. Figure 6 shows typical outcome. Physically, the 
R spectrum shows the frequency dependence of the correlation coefficient. The commonly used time-domain 
correlation data is unable to show this dependence. An examination of figure 6(b) shows that fluctuations 
around St = 0.6 remain well correlated (R = 0.9) after 0.5 D separation while those around St = 2.5 are less 
correlated (R = 0.3). Following is a discussion of the measurement uncertainty. The calculation of the 
coherence spectrum involves a division process which poses a problem with experimental data. The 
inevitable random noise present in data sets is greatly amplified in the final result. To reduce the impact, the 
auto- and cross-spectra were first smoothened by applying a Savitzky-Golay filter (ref. 1) before calculating 
the coherence spectrum. The uncertainty in R is relatively small in the well converged part of the cross- 
spectrum. As the separation between the probe volumes is increased the higher frequency part becomes 
progressively less correlated. This, at first, is reflected as an increase in the random error in the tail end of the 



Figure 6. — (a) Auto spectrum of turbulent density fluctuations at X(3.0D, 0.5D) and cross-spectrum 
from one probe at the same reference location and the other separated axially by 0.5D in M, = 1.4 jet. 
The dotted line show the smoothened spectra which was used to calculate coherence spectrum R 
shown in (b) using equation (10). 
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cross-spectrum. With increasing probe separation the high frequency part becomes submerged in the 


up at the high-frequency end. To avoid larger error the frequency range for which R<0.2 was not used for the 
subsequent analysis. 

Coherence spectra for increasing probe separation are shown in figure 7. For the most part such spectra 
decay with increasing probe separation. The rate of reduction is faster with increasing frequency, which 
reflects a shorter length scale for higher frequency eddies. The frequency range 0.15<St<0.8 in the shear layer 
data (fig. 7(a)), however; shows lasting correlation even when the probes are separated by four diameters. 
Frequencies above and below this range show quicker decay. The Kelvin-Hel mh oltz instability waves present 
at the jet shear layer are expected to be in the above frequency range. The increased correlation is reflective of 
longer length scales associated with these instability waves. At the centerline and beyond the end of the 
potential core the coherent structures are of lower frequency. This is reflected in R spectra shown in 


Coherence spectra are useful in calculating frequency dependent length scales. Towards this end R values 
for desired frequencies and for all probe separations were extracted and plotted in figures 8(a) and (c). Unlike 
the time-domain correlation r, the frequency domain coherence R does not become negative, therefore; 
problems associated with calculating integral length scale are avoided. Note that for certain frequencies 
coherence values become greater than unity. This is indicative of a growing level of fluctuation due to the 
turbulence production as discussed after equation 10. Nevertheless, R>1 does not affect integral length scale 
calculations. An exponential curve R = Ciexp(-C2 q x ) was fitted through data sets for each frequency; Ci and 
C 2 are fitted constants. To reduce the level of uncertainty data points close to the noise floor were left out 
of the curve fitting process. The integral length scale was determined as: 


An examination of data presented in this table provides many interesting observations, (a) As expected 
the integral length scale decreases with increasing Strouhal frequency. In other words, smaller eddies 
have higher frequencies, (b) The wideband length scale, traditionally used in noise source modeling, is 
representative of the mid-frequency range 0.7<St<1.0. (c) At the high frequency end, St>2, length scales 


Figure 7. — Variation of coherence spectra with increasing (indicated) probe separation. The reference 
probe is at (a) lip line X(3D, 0.5D), and (b) centerline X(8D, 0.D) of Mj = 0.95 jet. 


experimental noise floor. For these instances the coherence spectrum flattens out, and frequently starts to lift 


figure 7(b). 


00 


00 



( 22 ) 


The calculated values for the integral scales for different reference probe location are presented in table I. 


are similar in the shear layer and the centerline, this perhaps points towards a universal behavior of the 
high frequency eddies. However, at the low frequency end St<1.0 lip shear layer scales are longer than 
the centerline ones, (d) A rapid decrease in the length scale is observed radially across the lip shear layer. 


1.5 


(a) 


(b) 



o.o 


St = fD/Uj 20 


2.5 


0.5 


i.o St - fD/Uj 2.0 
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Figure 8. — Variation of coherence and cross-spectral phase with probe separation. Data 
from the same locations as in the earlier figure; Mj = 0.95, fixed probe at (a), (c) lip line 
X(3D, 0.5D), and (b), (d) centerline X(8D, 0.D). 

Returning back to figures 8(b) and (d), at a given frequency, the variation of the cross-spectral phase 
with spatial separation provides convective velocity (eq. 8) at that frequency. A linear fit through data 
points from different probe separations was used to reduce random error. The frequency dependent 
convective Mach numbers M c (St) = U c /Uj are shown in table II. An examination of this table shows the 
following, (a) The higher the Strouhal frequency the higher the convective Mach number. This is 
intuitively expected as the smallest eddies are expected to be convected close to the mean flow velocity, 
(b) M c is higher for fluctuations present in the centerline than at the lip line. This perhaps is a reflection of 
the higher mean velocity at the jet centerline, (c) M c decreases radially away from the jet centerline. 

The frequency and Mach number dependencies are further explored in figure 9. The length scales are 
normalized by Strouhal frequency L x St/D, yet the frequency dependence can not be totally eliminated. In 
general, this normalized length scales increases with frequency before attaining a near constant level. The 
same is true with the convective velocity. The shear layer data in figure 9(a) shows a low frequency range 
(St = 0.3 to 1.0) where length scales are unusually longer. The range varies with Mach number, 
nevertheless the longer length scales are indicative of Kelvin-Helmholtz instability waves which stay 
coherent over longer distances. The centerline data in figures 9(c) and (d) shows collapse of data points 
from all three jets under study. Note that reference points were selected such that they are approximately 
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St: fD/Uj St: fD/Uj 

Figure 9. — Variation of (a), (c) moving frame integral length scales and (b), (d) convective 
Mach number with Strouhal frequency in three different Mach number plumes. In (a) and 
(b) the reference probe was kept fixed at a shear layer position of X(3D, 0.5D); while in (c) 
and (d) it was moved to (8D, O.D), (10D, O.D), and (12D, O.D) in Mj = 0.95, 1.4, and 1.8 
jets respectively. 

1.2 times the potential core length away from the nozzle exit. The potential core length is a measure of the 
compressibility effect. The collapse of data shows that it should be used to scale the Mach number 
dependence of integral length scales. Note that the experimental uncertainty in L x at the highest resolved 
frequency is relatively higher than the lower frequency points. This is caused by a fast decay in coherence 
within a relatively small probe separation. Hence, only a few separation points can be used to estimate the 
high frequency length scales. 

Finally, figure 10 presents wave-number frequency spectrum from reference points at lip shear layers 
and centerlines of two different Mach number jets. For each reference point X, cross-correlation data 
were measured from n = 64 equally spaced, Aq x = 0.1D, probe separations (including the overlapped 
probe location) along the axial direction. The complex cross-spectrum G xr ,(X, q x , coi) was calculated for 
each probe separation following equation (20) (coi = 2nfi). Cross-spectral data for each frequency coi were 
extracted for all probe separations and a spatial Fourier transform was performed: 

Skffl( x > k x,m=®l) = y G xT 1 expfj^^-l k x m =— ^ — m = -^-,....+^ (23) 

~ 0 V n ) nAr| x 2 2 
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Figure 10. — Magnitude of the frequency-wave number spectrum measured from the indicated probe 
locations in the indicated Mach number plumes. The diagonal lines correspond to a convective velocity 
equal to the ambient sound speed. 


The subscript x in k x- IT1 denotes the axial wave number which is evaluated at m different points. The 
magnitude of each spectrum, normalized by (pj-p a ) 2 , p a is ambient density and Pj is the core density, is 
shown in figure 10. The staircase like appearance of the color plots is reflective of the coarse wave 
number resolution, Ak x = 1.56, achievable with the limited spatial data. Additionally, the minimum value 
in the color scale was selected to avoid plotting the noise floor. The utility of the k-eo spectrum to directly 
determine the radiative component of turbulent fluctuations was discussed earlier in the introduction. The 
part of the spectrum that lies to the left of the sonic speed line eo/k x = c® is capable of radiating at different 
observer angles. While most of the fluctuations in the subsonic Mach 0.95 jet are incapable of radiating to 
the far field, the supersonic Mach 1.8 jet shows a different scenario. With an average supersonic 
convective velocity, fluctuations in both the lip shear layer and the centerline of the jet are expected to 
radiate strongly along 0 = cos -1 (c<x/U c ) ~ 20° to 25° to the jet axis. The figure identifies that some parts of 
the low frequency fluctuations in both jets are capable of radiating at all angles. In this regard, compared 
to the shear layer, fluctuations from the centerline of the jet are found to possess larger radiative 
component of fluctuations. 


Summary and Conclusion 

Two point space-time correlations of density fluctuations were measured in subsonic and supersonic 
unheated plumes using a non-invasive Rayleigh scattering based optical diagnostics technique. The 
density fluctuations in the unheated plumes are solely due to the compressibility effect, and are expected 
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to be closely related to the noise sources. The two-laser beam Rayleigh scattering setup allowed point- 
measurements from two independently located probe positions in the plume. The system was designed 
and built from scratch using available laser, optical components, and electronics. Photo- electron counting 
electronics were used to measure the intensities of the scattered light collected from probe volumes. The 
count rates were directly related to the local air densities. To measure time history of fluctuating air 
density, counting was performed in contiguous time bins over long durations. Analysis of the count 
signals provided various statistics of space-time correlations. The particular goal of the present work was 
to determine the frequency dependent correlation scales in high velocity jets which were unavailable till 
now. Following is a summary of the primary observations: 

(a) Correlations were measured from the lip shear layer as well as from centerline and downstream of 
the potential core. The two sets show significant differences. At longer time delays significant negative 
correlations were measured in the shear layer data. This was mostly absent from the centerline data. The 
negative correlation prohibited realistic calculations of integral length scales in a fixed frame of reference. 
The time envelope of correlations from many different probe separations provided moving frame integral 
time scale and convective velocity. These values provide “wideband” or average property of all 
frequencies of turbulent fluctuations. 

(b) To obtain frequency dependent correlation, cross-correlation data were Fourier transformed and 
normalized by fluctuations measured at the overlapped probe location. The coherence spectrum provided 
frequency dependent correlation coefficients, while the phase of the cross-spectrum provided measures of 
convective velocities. As expected the convection velocity is found to increase with fluctuations 
frequency. The variation of the magnitude of coherence spectra with probe separation allowed for 
calculations of the frequency dependent length scales L(St) in a moving frame of reference. Such length 
scales were found to decrease progressively with frequency. 

A closer look into the shear layer data shows a range of frequencies, typically 0.1<St<1.0 (the range 
varied with Mach number and probe location) where integral length scales are longer than the neighboring 
frequencies. The longer length scales are believed to be associated with the Kelvin-Ffel mh oltz instability 
waves. 

(c) Data obtained from different Mach number jets were found to collapse when the reference probe 
locations were varied according to the length of the potential core. 

(d) By taking correlation data from a large number of spatial separations and by taking spatial Fourier 
transforms of cross-spectral data sets, the frequency-wave number eo-k spectrum was constructed. A 
comparison of such spectra from different Mach number jets showed a significant increase in the radiative 
part of turbulent fluctuations in supersonic jets. For all Mach number jets low frequency fluctuations from 
the centerline of the jet, compared to that in the shear layer, are found to possess larger radiative 
components. 
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